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ABSTRACT 

We report results from extensive new XMM- Newton observations of the ul- 
traluminous X-ray source (ULX) NGC 5408 X-l, one of the few ULXs to show 
quasi-periodic X-ray variability. We detect quasi-periodic oscillations (QPOs) in 
each of four new 100 ks each) pointings, expanding the range of frequencies 
and rms amplitudes observed from the source to 10-40 mHz and 10-45 %, respec- 
tively. However, similarly significant variations in the power-law photon spectral 
index, P, are not observed. We use the results of timing and energy spectral mod- 
eling to compare with the timing and spectral correlations seen in stellar-mass 
systems. We find that the qualitative nature of the timing and energy spectra of 
NGC 5408 X-l are very similar to stellar-mass black holes in the steep power- law 
state exhibiting Type-C QPOs. However, in order for this analogy to quantita- 
tively hold we must only be seeing the so-called saturated portion of the QPO 
frequency - photon index (or disk flux) relation. Assuming this to be the case, 
we place a lower limit on the mass of NGC 5408 X-l of > 800M©. Alternatively, 
the QPO centroid frequency is largely independent of the spectral parameters, in 
which case a close analogy of NGC 5408 X-l’s mHz QPOs with Type-C QPOs 
in stellar systems is problematic. Measurement of the source’s timing properties 
over a greater range of spectral parameters (in particular the spectral index) is 
needed in order to definitively resolve this ambiguity. We searched all the avail- 
able data for both a broad Fe emission line as well as high frequency QPO analogs 
(0.1-1 Hz), but detected neither. We place upper limits on the equivalent width 
of any Fe emission feature in the 6-7 keV band, and of the amplitude (rms) of 
a high frequency QPO analog of « 10 eV and ss 4%, respectively. 

Subject headings: Ultraluminous X-ray source : NGC 5408 X-l : Timing-spectral 
correlations : Quasi-periodic oscillations : Intermediate-mass black hole candi- 
date (IMBH) 


Astronomy Department, University of Maryland, College Park, MD 20742; dheeraj@astro.umd.edu 

2 Astrophysics Science Division, NASA’s Goddard Space Flight Center, Greenbelt, MD 20771; email: 
tod.strohmayer@nasa.gov 



1. Introduction & Background 


Ultraluminous X-ray sources (ULXs) are bright, point-like, off-nuclear extragalactic X- 
ray sources. Since their discovery over two decades ago (Fabbiano et al. 1989), the true 
nature of these sources has been a mystery. Their isotropic X-ray luminosities range from a 
few x 10 39 ergs s -1 to as high as 10 42 ergs s -1 . Luminosities at the high end of this range are 
found in a handful of sources which have been dubbed hyperluminous X-ray sources (HLXs, 
e.g., M82 X-l: Matsumoto et al. 2001, Kaaret et al. 2001; HLX ESO 243-49 X-l: Miniutti 
et al. 2006, Farrell et al. 2009). The apparent X-ray luminosities of these sources exceed 
the Eddington limit for any known stellar-mass black hole, hence the term “ultraluminous” . 

Based on their variability on short timescales (some ULXs are known to vary on timescales 
of the order of a few minutes) and high luminosities, some ULXs axe almost certainly pow- 
ered by accretion onto black holes (excluding the population of X-ray bright supernovae; 
Immler & Lewin 2003). The controversy, however, is with their mass range. Black holes, in 
principle, can exist in three mass ranges. Stellar-mass black holes (StBH: Mass ~ 3-50 M 0 ) 
are the remnants of stellar evolution of massive normal stars (stars excluding the POP III 
category). On the opposite end of the mass scale are the Super-massive black holes (SMBHs: 
Mass ~ 10 6 - 10 Mq) - which are currently theorized to form either via accretion onto “seed” 
black holes in the early universe (~ few Myrs after the Big Bang) with subsequent mergers 
(Volonteri et al. 2003; Davies et al. 2011) or via direct collapse of dense gas (Begelman et 
al. 2006). Their exact formation scenario is still somewhat uncertain. Finally, one has the 
Intermediate-mass black holes (IMBHs: Mass range 100 - few x 1000 Mq), whose formation 
mechanisms are the subject of considerable debate. There are a number of processes whereby 
an IMBH can form. For example, they can be formed at the centers of globular clusters via 
accretion onto a massive black hole (Mass ~ 50 Mq) over a time period of 10 10 years (Miller 
et al. 2002), by accretion onto a primordial black hole (e.g., Mack et al. 2007), or by mergers 
of seed black holes. 

The off-centered nature of ULXs rules them out as candidates for SMBHs. Given the 
ages of their host galaxies, if they were SMBHs, they should have sunk to the centers of 
their galaxies due to dynamical friction (Miller & Colbert 2004). Clearly, this is not the 
case. This leaves the possibility of their being either StBHs or IMBHs. Theories supporting 
the idea of a StBH as the central engine often employ either beamed X-ray emission or 
some sort of super-Eddington emission mechanism to explain the observed high luminosities. 
Current ideas include geometric beaming (King et al. 2001), relativistic beaming similar 
to the phenomenon of blazars (Hording et al. 2002), and super-Eddington emission as 
a result of a photon-bubble instability in radiation-dominated accretion disks (Begelman 
et al. 2002). More recently, based on high-quality XMM-Newton spectra of 12 ULXs, 
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Gladstone et al. (2009) have suggested that ULXs might be accreting StBH systems in a new 
ultraluminous accretion state (an addition to the current classification of black hole accretion 
states: thermal (high/soft), hard (low/hard) Sz steep power law (very high) as in McClintock 
et al. 2006), where accretion is again super-Eddington. Alternatively, it has been suggested 
that these sources might be powered by gravitational energy release from accretion onto 
an IMBH (Colbert & Mushotzky 1999). This presents a relatively simpler solution, where 
the energetics proceed similarly to that of the StBH systems (at sub-Eddington accretion 
rates) with the corresponding physical properties scaling appropriately with the mass of the 
compact object. It is also possible that ULXs constitute an inhomogeneous population of 
both StBHs and IMBHs. 

It is now known that certain characteristic timescales in accreting black hole systems 
scale with the mass of the black hole. For example, McHardy et al. (2006) have established 
that in soft-state StBHs and AGN, the break timescale of the power density spectrum (PDS), 
the mass of the black hole and the luminosity are strongly correlated. Their work was 
further extended to the hard-state StBHs by Hording et al. (2007). These studies strongly 
demonstrate that certain physical timescales of accreting compact objects scale directly with 
the mass of the source (after accounting for differences in accretion rates, i.e., luminosity). 
Therefore, timing studies can play a key role in extracting valuable information about the 
mass of the compact source, which is crucial in the case of ULXs in order to better understand 
the nature of the emission processes that result in such high apparent luminosities. Of 
particular interest here are the low frequency quasi-periodic oscillations (LFQPOs). Within 
the context of StBHs, the LFQPOs (those in the range ~ 0.1-15 Hz) are broadly classified 
into three categories based on their properties and the overall nature of their PDS (Casella 
et al. 2005). The PDS with type- A and B QPOs is characterized by weak red noise (noise 
at the low frequency end of the PDS) with the type-A QPOs occurring with relatively low 
coherence (quality factor, Q: centroid frequency /FWHM, < 3) compared to the type-B (Q 
> 6). Finally, one has the type-G QPOs which are most relevant to the present work. The 
PDS accompanying these QPOs in StBH systems can be described by a flat-topped, band- 
limited noise breaking to a power law with the QPOs evident on the power-law portion of 
the spectrum, close to the break. 

Further, the type-C QPOs are fairly coherent with the quality factor, Q from 5-15 and 
amplitudes (% rms) ranging from 2-20. In StBHs, they are known to occur in the frequency 
range from f» 0.1-15 Hz. In light of the work unifying the stellar-mass and the super-massive 
black holes, if some ULXs were to host IMBHs then the qualitative behavior of their PDS 
should be comparable to StBHs with the characteristic variability times scaling according 
to the mass of the putative IMBH. This idea has been explored by a number of authors to 
search for the "QPO analogs” in ULXs; and QPOs have now been detected in a handful of 
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them. The QPOs detected in M82 X-l (Strohmayer h Mushotzky 2003), NGC 5408 X-l 
(Strohmayer et al. 2007) and NGC 6946 X-l (Rao et al. 2010) resemble the type-C QPOs, 
while those detected in the M82 source X42.3+59 more closely resemble the type-A or -R 
QPOs (Feng et al. 2010). The crucial difference here is that the QPO centroid frequencies 
of the ULX sources appear to be scaled down by a factor of a fewx(10 - 100) (~ few mHz) 
compared to the LFQPOs in StBHs. 

Under the assumption that the ULX mHz QPOs are analogs of the LFQPOs in StBHs, it 
is reasonable to assume that their characteristic timescales/frequencies (e.g., QPO centroid 
frequencies, PDS break frequencies) scale with the mass of the accreting source and vice 
versa. However, the type-C QPOs occur with a wide range of centroid frequencies (0.1 — 15 
Hz) in StBH systems. Therefore, timing information alone is not sufficient to accurately 
estimate ULX masses in this way, but combining timing and spectral information has proven 
to be a valuable tool. For example, in StBH systems, the power-law photon index and disk 
flux are correlated with the QPO centroid frequency. The general trend is an increase in 
power-law photon index and disk flux with the QPO centroid frequency (e.g., Sobczak et al. 
2000; Vignarca et al. 2003; Shaposhnikov & Titarchuk 2009), with evidence for saturation 
(constancy of the power-law photon index and disk flux with a further increase in QPO 
centroid frequency) beyond a certain frequency. 

Using a reference StBH system with a measured QPO centroid frequency - photon- 
index relation, one can then scale the QPO centroid frequencies detected in LTLX systems at 
a given power-law spectral index to get an estimate of the mass. For example, with archival 
X M M-Newton data from M82 X-l, Dewangan et al. (2006) extracted the energy and power 
spectra of the source. The PDS was strikingly reminiscent of a StBH with type-C QPOs, i.e., 
flat-top noise breaking to a power-law with QPOs on the power-law portion of the spectrum. 
However, the respective timescales were scaled down by a factor of ~ 10. The QPO in this 
case was centered around a frequency of ~ 114 mHz and the energy spectrum had a power 
law photon index of ~ 2.0. Using the QPO frequency - photon index correlations from two 
StBH reference sources, GRS 1915+105 and XTE J1550-564, they estimated the mass of the 
ULX in M82 X-l by scaling its QPO centroid frequency (~ 114 mHz) at the given photon 
index (~ 2.0). They estimated the mass of the black hole to be in the range 25-520A/©. 
Similar scaling arguments were used by Rao et al. (2010) to estimate the mass of the black 
hole in the ULX NGG 6946 X-l to be in the range (1-4) x 1000 M©. Based on both the 
PDS and the energy spectrum of NGC 5408 X-l, Strohmayer et al. (2009) (S09) argued 
that the source behavior was consistent with the steep power-law state (SPL) often seen 
in StBHs. They compared the available data from NGC 5408 X-l to five different StBH 
reference sources and estimated the mass of the black hole to be a few x 1000 M©. Feng 
et al. (2010) detected 3-4 mHz QPOs from the ULX X42.3+59 in M82 and identified them 
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as either type A/B analogs of StBHs. They estimated the mass of the black hole to be in 
the range 12,000-43.000 M© by scaling the QPO frequency to that of the type A/B QPOs 
in StBHs. 

It is important to note that the above scaling arguments have several caveats. First, 
the mass estimates were established under the assumption that the mHz QPOs detected 
from these ULX systems are analogous to a particular type of LFQPO detected in StBH 
systems, i.e., A, B or C. These identifications were based on the qualitative nature of the 
power spectrum (PDS) alone in some cases and both the PDS and the energy spectrum in 
the case of NGC 5408 X-l. Second, in the case of the ULX mHz QPOs, the observed range 
of QPO centroid frequencies and photon indices has been limited. To gain a more secure 
identification one would like to see the QPO frequencies and photon spectral indices correlate 
in a similar fashion as for the StBHs. 

In this paper we further explore these issues using extensive new observations of NGC 
5408 X-l with XMM-Newton. We describe the properties of the mHz QPOs over a wider 
range of centroid frequencies than previous data allowed. More specifically, we study the 
nature of the mHz QPOs from NGC 5408 X-l through a systematic search for timing - 
spectral correlations similar to those seen in the StBHs. The article is arranged as follows. 
In §2 we describe the data used for the present study. In §3 we show the results from our 
timing analysis, while in §4 we give details of the energy spectral analysis. In §5 we describe 
our search for timing - spectral correlations similar to those seen in accreting StBH systems. 
The search is conducted using two different spectral models, a phenomenological model of a 
multi-colored disk + power-law and then with a model describing Comptonization by bulk 
motion (Titarchuk et al. 1997) that has been used previously to derive black hole masses 
from QPO scaling arguments (Shaposhnikov & Titarchuk 2009). Finally, in §6 we discuss 
the implications of our results on the mass of the black hole in NGC 5408 X-l. 


2. XMM-Newton Observations 

XMM-Newton has now observed NGC 5408 X-l on multiple occasions. We use data 
from six of the most recent observations (~ 100 ks each) spread over a time span of five 
years (2006-2011). The details of the observations are outlined in Table 1. Results from 
the first two observations (in 2006 and 2008, respectively) were summarized in Strohmaver 
et al. (2007) (S07) k Strohmayer et al. (2009) (S09). Using the 2006 data, S07 reported 
the first detection of quasi-periodic variability from this source. S09’s analysis of the 2008 
observation again showed evidence for the presence of quasi-periodic oscillations. However, 
the most prominent QPO during the 2008 observation was at a lower frequency (QPO at 
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~ 10 mHz compared to 20 mHz in 2006). Further, S09 noted that the disk contribution to 
the total flux and the power-law index of the energy spectrum decreased slightly compared 
to its state in 2006. This is analogous to a trend often seen in StBH systems where the 
disk flux and the power law index of the energy spectrum positively correlate with the 
centroid frequency of the most dominant QPO (e.g., Viagnarca et al. 2003, Sobczak et 
al. 2000). These findings were used to propose additional observations with XMM-Newton, 
with the goal of detecting QPOs over a range of frequencies and hence to further explore for 
correlations between timing and energy spectral properties. A large program was approved 
for Cycle 9 (PI: Strohmayer), and four observations (~ 100 ks each) were made under this 
program. Two of the observations were carried out in 2010 (2010a & 2010b) and the rest 
in 2011 (2011a &; 2011b). Here we present results from these new observations as well as 
a reanalysis of the earlier pointings, so as to facilitate a consistent comparison of all the 
available data. 

For the present work, we use the data acquired by the European photon imaging camera 
(EPIC), i.e, both the pn and MOS to get a higher signal-to-noise ratio in the power and energy 
spectra. We used the standard SAS version 11.0.0 to reduce the images and filtered event- 
lists from all the EPIC data. The standard filter with (PATTERN<=4)> f° include only 
single and double pixel events, was applied to the event-lists and events only in the energy 
range 0.3-12.0 keV were considered for further analysis. Background flaring was prominent 
for brief periods during certain observations. The power and energy spectra were carefully 
extracted taking into account both the background flaring and the instrumental dead time 
effects (specific details in the next sections). In all the observations the source was easily 
detectable and we did not face any source confusion problems. We extracted source events 
from a region of radius 32” centered around the source. This particular value was chosen 
to roughly include 90% of the light from the source (estimated from the fractional encircled 
energy of the EPIC instruments). A background region, free of other sources, was extracted 
in a nearby region. The size of the background region was chosen to be consistent with the 
source region. Further, the size of the source and the background region was chosen to be 
consistent for all the six observations (i.e., 32”). We present the specific details of the timing 
and spectral analysis in the following sections. 


3. Timing Analysis 

We first produced source and background light curves to assess the quality of the 
data. Shown in the left panels of Figure 1A and IB are the background count rates 
(pn+ M OS 1 + M OS 2) for each of the six observations. Overlaid are the time intervals during 
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which a given instrument (pn, M OS 1/2 ) was continuously active for more than 5 ks, i.e., 
all the instrumental good time intervals (GTIs) longer than 5 ks. For a given instrument, a 
horizontal line (offset to an arbitrary value for each instrument) indicates the active time, 
while a vertical line marks the beginning or the end of a continuous data interval. Such an 
insight is important as we are using both the pn and the MOS data to achieve higher count 
rates compared to, say, the pn alone. In other words, this ensures that the combined pn and 
MOS event-lists contain only events corresponding to the times during which all the three 
instruments were active. Furthermore, inspection of the light curves reveals background 
flaring at either the beginning or the end, or both, of all the observations. Therefore, the 
final combined (pn+MOSl+MOS2) light curves were extracted by taking into account the 
GTIs of all the instruments and excluding the periods of background flaring. The filtered 
light curves (100s bins) of the source and the non-flaring background in the energy range 0.3 
- 12.0 keV are shown in the right panels of Figures 1A and IB. The combined mean count 
rates (0.3-12.0 keV) during each observation are listed in Table 1. Notice that the gaps in 
the data are either due to exclusion of background flaring times, dead time intervals or the 
absence of GTIs longer than 5 ks. 

It is evident from the light curves (right panels of Figure 1A and IB) that the source 
varies significantly during all the observations. To quantify the variability, we construct a 
power density spectrum from each of the light curves. To achieve higher signal-to-noise in 
the PDS, we break the GTIs into shorter segments (S/N cc -/Number of individual spectra, 
e.g., Van der klis 1989) and derive an average power spectrum. For a given observation, 
the size and the number of such segments are shown in the first row of Table 2. Figures 
2 A and 2B show the PDS and the best-fitting model (thick solid line) for each of the six 
observations in two different energy bands. Shown in the left panel is a PDS in an energy 
band in which a QPO is detected with high statistical significance and in the right panel is 
a PDS derived from photons in the energy band of 1.0-10.0 keV (the reason for constructing 
power spectra in two different energy bands is discussed in the next paragraph). All the 
power spectra shown here use the so-called Leahy normalization, with the Poisson noise 
level being 2 (Leahy et al. 1983). Clearly, in each spectrum the power rises below ~ 0.1 Hz 
with evidence for a QPO in the range of 10-40 mHz and essentially Poisson noise at higher 
frequencies. 

To quantify this behavior, we fit a bending power law to the continuum and a Lorentzian 
to model the QPO (Belloni et al. 2002). The model fits well with the reduced y 2 ~ 1 (last 
row of Tables 2 & 3) in three cases (2008. 20 11 A, 201 IB) and gives acceptable fits with the 
reduced y 2 in the range 1.3-1. 5 (last row of Tables 2 & 3) in the other three cases (2006, 
2010A, 2010B). A careful analysis of the residuals in the latter three cases indicates that 
multiple weak features contribute significantly to increase the overall y 2 . For example in 
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the case of the PDS derived from the 2010A data (bottom two panels of Figure 2A), the 
weak QPO-like feature at ~ 55 mHz and the excess at ~ 0.2 Hz contribute about 30 to the 
total x 2 - However, their individual statistical significance is rather low. Nevertheless, for 
the purposes of analyzing the QPO properties and studying the overall qualitative nature of 
the power spectra, the fits are entirely adequate. The best-fitting model parameters (derived 
from a fit in the frequency range 1.0 mHz - 0.5 Hz) for each of the observation are shown in 
Table 2 & 3 for the two different energy bands (highlighted in the tables). Also shown are the 
X 2 /degrees of freedom (dof) values for each of the fits along with the x 2 /dof corresponding 
to the continuum model (in braces). The change in x 2 serves as an indicator of the statistical 
significance of the QPOs. For a given observation, we choose two different energy 
bands: the first energy band (second row of Table 2) is the bandpass in which 
the QPO is detected with a very high confidence and the second bandpass of 1.0- 
10.0 keV was chosen to consistently compare the properties of the QPO across 
all the observations (overall variability is energy dependent). We confirm the overall 
qualitative nature of the PDS, that is, flat-topped noise breaking to a power law with a QPO 
on the power law portion (near the break). This is consistent with StBHs showing type-C 
QPOs but with the characteristic frequencies scaled down by a factor of a few x 10. This is in 
agreement with the results reported by S07 and S09. However, we cannot comment on 
the possibility of the QPO pairs as all our PSD (Figure 2A & 2B) are averaged 
to the extent of smearing out weak features (to clearly identify the strongest 
QPO). 

One of the main goals of the present work is to better characterize the mHz QPOs 
seen from NGC 5408 X-l within the context of the known classes of LFQPO seen in StBHs 
(Casella et al. 2005). An important diagnostic for this purpose is the study of the variation 
of the rms amplitudes of these QPOs with their centroid frequency. In the case of LFQPOs 
from StBH systems this behavior is fairly well-established (e.g., Revnivtsev et al. 2000; 
Sobczak et al. 2000; Vignarca et al. 2003; McClintock et al. 2006). The typical behavior 
of the StBH LFQPOs, in the energy range 2.0-20.0 keV (the nominal RXTE bandpass), 
can be described as follows: as the source evolves along the low/'hard state towards the 
SPL state, the QPO frequency increases and the rms amplitude increases. Before reaching 
the SPL state, the source traverses an "intermediate” state in which the QPO frequency 
continues to increase, but the rms amplitude decreases. Upon reaching the SPL state the 
correlation tends to break down, showing more scatter in the rms amplitude (see Figure 11 
in McClintock et al. 2009). 

The X-ray variability of NGC 5408 X-l is known to depend on energy (Middleton et al. 
2011). The overall rms variability increases with an increase in the energy of the photons at 
least up to 2.5 keV (Given the poor signal-to-noise at higher energies, it is not clear whether 
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the variability strength levels off or decreases). This is similar to the energy dependence of 
LFQPO (« few Hz) detected from the galactic micro-quasar GRS 1915+105 (Rodriguez et 
al. 2004), where the QPO amplitude increases with energy, before rolling over. We used the 
QPO fit parameters from all the observations to explore the dependence of the QPO rms 
amplitude with centroid frequency. To be consistent across all the observations, we chose 
an optimum bandpass of 1.0-10.0 keV to ensure good statistics and significant detection of 
the variability (the variability of the source increases with higher energy while the 
countrate decreases). In addition, to explore whether our results might depend on energy, 
we derived the QPO properties in a different energy band (though with some overlap of the 
optimum band defined above) of 0.5-5. 0 keV. The rms amplitude versus QPO frequency 
results for the 1.0-10.0 keV and 0. 5-5.0 keV are shown in the left and the right panels of 
Figure 3, respectively (compare with Figure 11 of McClintock et al. 2009). Given that the 
QPO amplitude is higher at higher energies (Middleton et al. 2011), not all power spectra 
derived from the lower energy band (0.5 - 5.0 keV) yielded statistically significant QPO 
features. This is why the right panel of Figure 3 has one less data point. Clearly, the 
variation is similar in both the energy bands. We further compare with the results from 
StBHs in §5. 

Another prime driver for the long observations of NGC 5408 X-l was to search for high 
frequency QPO (HFQPO) analogs of StBHs. These QPOs are observed in the range from 
about 50 - 450 Hz in StBHs (McClintock & Remillard 2006; Strohmayer 2001a). Simple mass 
scaling arguments would suggest that in IMBHs they would be expected in the frequency 
range of 0.1-1 Hz, with a low rms amplitude of ~ 2%. We do not detect any obvious QPO- 
like feature at these frequencies in any of the individual observations. Further, given that 
HFQPOs appear to be reasonably stable in frequency in StBHs (e.g., Strohmayer 2001b), we 
averaged the power spectra derived from the individual observations to improve the signal- 
to-noise, however, this did not lead to a detection. We estimate the upper limit on the rms 
amplitude of a QPO-like feature in the frecjuency range of 0. 1-1.0 Hz to be ~ 4%. 


4. Energy spectral Analysis 

We fit the X-ray spectra of NGC 5408 X-l using the XSPEC (Arnaud 1996) spectral 
fitting package and the EPIC response files were generated using the arfgen and rmfgen 
tools which are part of the XMM-Newton Science Analysis System (SAS) software. Since, a 
primary goal of the present work is to search for timing - spectral correlations similar to those 
seen in accreting StBHs, we began our analysis by characterizing the energy spectra using 
the same phenomenological models often used to describe the X-ray spectra of accreting 
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StBHs. In terms of XSPEC models, we used diskpn + power-law. We also explored the hmc 
model, Comptonization by bulk motion (Titarchuk et al. 1997), because it has been used 
to derive mass estimates from QPO frequency scaling arguments, and reference spectral - 
timing correlations have been derived for a significant sample of StBHs using the parameters 
derived from this model (see Shaposhnikov & Titarchuk 2009). We start by describing some 
of the specifics of our data extraction and reduction methods. 

Since we are interested in exploring the correlations between the timing and the spectral 
behavior, we elected to extract energy spectra from the same time intervals as those used to 
extract the PDS (see the first row of Tables 1 & 2). Such synchronous measurements will 
tend to minimize any offsets that could be induced by variations in the system properties 
(both spectral and timing) within a given observation. In addition to the standard filters 
described in §2, the (FLAG=—0) filter was imposed to get the highest quality spectra. The 
combined ( pn+MOS ) average count rate of the source varied in the range 1.19 - 1.46 cts 
s -1 over all the observations. Given such relatively low count rates, source pileup was not 
an issue in any of the six observations. Since we are using both the pn and the MOS data, 
we individually reduce the energy spectra from each of the detectors using the same good 
time intervals. This gave us three energy spectra (pn, MOS1, MOSS) for each of the six 
observations. A given model was fit simultaneously to all the three spectra to derive the 
tightest constraints on the best-fit parameters. Both the pn and the MOS spectra were 
re-binned to ensure a minimum of 25 counts in each bin. Given the high number of total 
counts (ss (6-12) x 10,000), this gave us high-quality spectra in each case. 

We began by fitting the spectra with a multi-colored disk ( diskpn in XSPEC, Gierlinski 
et al. 1999) + power law ( powerlaw in XSPEC) + an X-ray hot plasma ( apec in XSPEC) as 
in S09. We used the Tuebingen-Boulder ISM absorption model ( tbabs in XSPEC) to account 
for interstellar absorption. The hydrogen column density was set as a free parameter. This 
model is very similar to that often used to model the spectra of StBHs. Indeed, it is common 
to use this or similar phenomenological models to determine the spectral state of a StBH. 
However, here in addition to a thermal disk and a corona, there is an X-ray emitting plasma 
in the model. This was realized by S07 after identifying systematic emission-like features in 
the residuals of the diskpn + power-law fit to the 2006 data. We confirm their result that 
the x 2 improves significantly (by ~ 40-100 over all the six observations) with addition of the 
apec component to the standard multi-colored disk + power-law model. Such a feature is not 
exclusive to the ULX in NGC 5408, as NGC 7424 (Soria et al. 2006c) and Holmberg II X-l 
(Dewangan et al. 2004) also show evidence for the presence of an X-ray emitting plasma. 
In addition, the high-resolution VLA radio observations at 4.9 GHz of the counterpart of 
the X-ray source in NGC 5408 X-l revealed a radio nebula of « 40 pc extent (Lang et al. 
2007). Based on the value of the power-law index of the radio spectrum, it was suggested 
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that the radio emission is likely from an optically thin synchrotron emitting gas. Given this, 
it is possible that the putative X-ray hot plasma is coincident with the radio nebula. The 
best-fitting parameters of the above model for each of the six observations are shown in 
the top panel of Table 4. We note that the X-ray temperature of the plasma has 
remained constant across the suite of our XMM-Newton observations. 


We also fit the spectra with a model based on Coinponization by bulk motion (Laurent 
& Titarchuk 1999, bmc in XSPEC). As mentioned above, this model has been used within the 
context of quantifying timing - spectral correlations in StBHs (Shaposhnikov & Titarchuk 
2009). The underlying physical basis of this model is similar to that of the multi-colored disk 
+ power-law model, i.e. , the presence of a thermal disk and a hot electron corona. However, 
in this case the power-law is produced by the Comptonization of soft photons (from the disk) 
within a converging inflow onto the black hole. In this model, the normalization parameter 
derived from the fit mimics the disk flux from the diskpn component in the other model we 
have employed. We also include the apec component in this case. The best-fitting model 
parameters including the absorbing hydrogen column density using the tabs in XSPEC for 
each of the six datasets are shown in the bottom panel of Table 4. 


It can be noted straight away that, for a given observation, the plasma temperatures 
{apec) found from fitting the two spectral models are consistent with each other (within the 
error bars). The reduced x 2 values for the diskpn+pouilaw+apec and the bmc+apec model are 
in the range of 1.02-1.24 and 0.99-1.15 (~ 600-900 degrees of freedom: last column of Table 
4), respectively. Clearly, they both give acceptable values of reduced y 2 . The latter model 
( bmc+apec ) fits the data slightly better (see Table 4) than the simple diskpn+apec+powlaw 
model. However, the improvement in y 2 is not statistically significant in all the cases and 
therefore it is not possible to rule out one model over the other based on the reduced y 2 

cu r va t- m ?& ) - H h gh t- better describe- the high - c[ua lit y ^ d a t ft- e f-HL X s' (^g-.- r Gl adston e ct al. 2009). 
A sample energy spectrum (using the 2010B data) is shown in Figure 4. Also shown in the 
figure are the ratios of the data to the folded model in the two cases {diskpn+apec+powlaw 
and bmc+apec). We use the spectral parameters derived here and the timing properties from 
the previous section to search for correlations between the two. The results of w'hich are 
presented in the next section. 


The measurement of fluorescent lines from elements like Fe, Ni, Cr, Ca etc. is another 
important probe of the vicinity of the X-ray emitting region. In principle, these reflection 
features can be produced by hard X-rays, presumably from the Comptonizing corona, irra- 
diating cold material, e.g.. in the disk. Among all such emission lines with their rest-frame 
energies in the XMM-Newton bandpass (0.2-10.0 keV), the iron feature is the most promi- 
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nent (e.g., Matt et al. 1997). Given the subtle changes in the energy spectra over the suite 
of observations (see Table 4), it does not seem unreasonable to combine all the data while 
searching for weak emission features. We searched for iron emission line features in the 
energy range of 6. 0-7.0 keV, but did not detect any obvious features. Using an effective ex- 
posure of « 540 ks (all datasets, see Table 1), we were able to place tight constraints on the 
equivalent width of any feature in the energy range of 6. 0-7.0 keV. Assuming an unresolved 
narrow feature of width 10 eV, the upper limit (90% confidence) on the equivalent width of 
an emission line at 6.4 keV is 5.4 eV. We then check for the possibility of broad emission 
lines in the energy range 6. 0-7.0 keV. The upper limits on the equivalent width assuming a 
broad emission line of width 0.3 keV and 1.0 keV are 11.6 eV and 11.0 eV, respectively. 

The apparent weakness of iron emission from the disk can be due to a number of factors: 
an accretion disk that is completely ionized, scattering of the reflected iron-line photons in 
an optically thick corona (Matt et al. 1997), especially when the corona is a thin layer above 
the accretion disk or a low iron abundance in the accretion disk (Matt et al. 1997). Further, 
a high inclination angle for the accretion disk combined with one or more of the above factors 
can decrease the intensity of the iron-line. Recently, Gladstone et al. (2009) have analyzed 
energy spectra of a sample of ULXs (including NGC 5408 X-l) and have suggested that 
these sources might be operating in a new accretion state, with rollover at higher energies 
(> 2.0 keV) as a characteristic signature of such a state. Thus, the nature of the spectrum 
might be such that there are not enough photons at higher energies to generate a detectable 
iron-line fluorescent feature (Cole Miller, private communication). 


5. Timing - Spectral Correlations 

A primary goal of the new observations of NGC 5408 X-l is to further test the pre- 
liminary classification of the mHz QPOs detected from NGC 5408 X-l as the analogs of 
the 0.1-15 Hz Low-Frequency, type-C QPOs (LFQPOs) detected from accreting StBHs. If a 
strong connection can be demonstrated, then scaling of the characteristic frequencies (e.g., 
QPO centroid frequency) can be better justified to estimate the mass of the black hole within 
this source. 

As mentioned earlier, scaling relations using the PDS break frequency have been suc- 
cessfully tested using power spectra of both stellar and super-massive black holes of known 
mass (McHardy et al. 2006 and Hording et al. 2007). Based on the qualitative nature of 
the PDS and the energy spectra derived using the 2006 and 2008 data, S09 suggested that 
the mHz QPOs from NGC 5408 X-l may indeed be LFQPO analogs. More specifically, 
they proposed that its X-ray state was analogous to that of a StBH in the so-called steep 



power-law (SPL) state exhibiting type-C QPOs, but at the same time emitting a few x 10 
higher X-ray flux. Having analyzed the data from all six observations (§3 k 4) we confirm 
that the derived power- and energy-spectral parameters are qualitatively consistent with the 
source being in a SPL state exhibiting type-C QPOs (based on the state descriptions in 
McClintock k Remillard 2006). A summary of the working definition of the SPL state with 
type-C QPOs and the PDS k energy spectral parameters of NGC 5408 X-l are given in 
Table 5. 

A further step in understanding the nature of these QPOs is to determine if they show 
the same evolutionary behavior in the QPO parameters as those exhibited by the StBH 
LFQPOs. One such characteristic behavior is the QPO rms amplitude versus centroid fre- 
quency relationship. Our results for this relation from NGC 5408 X-l are shown in the left 
and right panels of Figure 3 in two different energy bands (1.0-10.0 keV and 0. 5-5.0 keV, 
respectively). Comparing these plots with those from StBHs, as in Figure 11 of McClintock 
et al. (2009), the results appear to be at least qualitatively consistent with the behavior of 
SPL data from some StBHs (see for example, the green triangles in Figure 11 of McClintock 
et al. 2009). However, because a tight correlation appears to break down in the SPL state 
it is hard to conclude definitively that the NGC 5408 X-l behavior is exactly analogous. 

Another important signature of the Type-C LFQPOs in StBHs is their dependence on 
the energy spectral features, i.e., the timing - spectral correlations. Using the X-ray state 
classifications as in McClintock et al. (2006), LFQPOs in StBHs are usually detected in 
the low/hard state, SPL state and during the transition between these two states. Within 
the context of the multi-colored disk + power-law model parameters, it is known that the 
centroid frequency of these LFQPOs is correlated with the disk flux and the photon index of 
the power-law component. In StBHs, the typical behavior is that the QPO centroid frequency 
is positively correlated with the disk flux and the photon index of the power law. At a certain 
higher QPO centroid frequency (« 5-10 Hz) the relationship seems to flatten (saturate) and 
even appears to start reversing in some cases (Vignarca et al. 2003). Furthermore, the 
correlation is not exactly the same for all StBHs, i.e., the value of the slope of the correlation 
and the photon index/disk flux at saturation are different for different sources and can be 
different for the same source in a different outburst (Vignarca et al. 2003; Shaposhnikov & 
Titarchuk 2009). In Figure 5 we show plots of several derived spectral parameters and fluxes 
(using the diskpn + power-law fits) as a function of the QPO centroid frequency. These can 
be compared with similar plots derived from StBH systems (see, for example, Figures 8 k 9 
in McClintock et al. 2009; Figure 10 of Vignarca et al. 2003). Again, there appears to be 
qualitative consistency between the behavior exhibited by NGC 5408 X-l and StBHs, but the 
full range of behavior is not yet seen, again because the range of spectral variations in NGC 
5408 X-l (most notably in the power-law index) is too modest. Additional measurements of 
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QPO properties at lower values of the power-law photon index, for example, could provide 
a more definitive test. 

One can also investigate the spectral - temporal correlations in the context of the spectral 
parameters derived using the bmc model. Indeed, Shaposhnikov Sz Titarchuk (2009) explored 
such correlations in detail for a sample of StBHs using RXTE data. They also developed a 
set of fitting functions to quantify the observed correlations and used these to obtain mass 
estimates by scaling arguments. Figure 5B compares the results from NGC 5408 X-l with 
those from StBHs using the bmc spectral fits. We used the published data from Shaposhnikov 
& Titarchuk (2009) to produce representative plots for a StBH (XTE J1550-564, top panels). 
Shown from left to right are the QPO centroid frequency versus the power-law photon index 
(F), the bmc normalization (disk flux) versus F and the Comptonized fraction, / versus F, 
respectively. The data points are color coded to highlight the low/hard state (blue) and the 
SPL state (green). We used the data from the 1998 outburst of XTE J1550-564 which clearly 
demonstrates the typical behavior of StBHs (top panels). Comparisons of these plots with 
those in, for example, Shaposhnikov & Titarchuk (2009) lead to similar conclusions as with 
the results in Figure 5A. There is qualitative consistency with the Type-C identification, but 
the classification is not definitive. In order to make a more secure association of the QPOs 
in NGC 5408 X-l with Type-C (SPL state) QPOs in the StBHs, we need, at a minimum, to 
obtain observations of the source over a greater range of power-law spectral indices. 


6. Discussion 

It has been established that black hole masses scale with the break frequency of their 
PDS (McHardy et al. 2006, Hording et al. 2007) and this relation is known to hold over 
6 orders of magnitude in mass, i.e., from StBH to super-massive black holes. Furthermore, 
at least in StBHs, the break frequency of the PDS is known to strongly correlate with the 
centroid frequency of the LFQPOs (Wijnands et al. 1999). Therefore, it is reasonable to 
assume that the mass of the black hole scales with the centroid frequency of the LFQPOs. 
This can be especially useful in cases where black hole mass measurements are otherwise 
difficult, viz., in ULXs. However, an important caveat is to ensure that the QPOs used for 
comparison are similar in nature. 

Our new results confirm the earlier result from S09 that the qualitative nature of the 
PDS and the energy spectrum of NGC 5408 X-l are very similar to that of StBHs in the 
SPL state (see Table 5). At the same time, the characteristic timescales within the PDS are 
lower by a factor of ~ 100, while the X-ray luminosity is higher by a factor of a few x 10, 
when compared to a typical StBH in such a state. However, based on the scaling 



between the break frequency of the PSD and the centroid frequency of the 
LFQPOs (Centroid frequency of LFQPOs « 12xbreak frequency of the PSD), 
the association of the mHz QPOs from NGC 5408 X-l with the LFQPOs from 
StBHs has been questioned (Middleton et al. 2011) 

A striking feature of the behavior of NGC 5408 X-l is that the energy spectrum of the 
source has remained roughly constant over the current suite of observations, while the prop- 
erties of the most prominent QPO (e.g., rms amplitude, centroid frequency) have changed 
significantly. These results appear consistent with the source being in a spectral state sim- 
ilar to the SPL state in StBHs. More specifically, in StBHs, certain spectral properties like 
the disk flux, the photon index of the power-law correlate positively with the QPO centroid 
frequency up to a certain frequency (fs 2-10 Hz) beyond which the relationship turns around 
or remains roughly constant. Detailed analysis of the timing - spectral correlations in StBH 
sources such as XTE J1550-564 and H1743-322 has revealed a comprehensive picture (Mc- 
Clintock et al. 2009) of their behavior. Tracking these two sources as they evolved from the 
low/harcl state into the SPL state, it was realized that the spectral properties like the disk 
flux and the power- law photon index are tightly correlated with the QPO centroid frequency 
during the low/hard state and during a transition phase intermediate between the low/hard 
and the SPL state; And the relationships seem to saturate/breakdown (spectral properties 
remain roughly constant with increase in the QPO frequency) as the sources entered the 
SPL state. It may be that we are in fact seeing NGC 5408 X-l in a similar state, that is, we 
may be seeing the “saturated” portion of the relation. 

Given that the overall accretion timescales change with the mass of the black hole 
(timescales increase as the mass increases from StBHs to super-massive black holes, e.g., 
McHardy et al. 2006), it is possible that we are seeing a similar timing - spectral behavior, 
and that the apparent longevity of the source in the SPL state might also be due to its 
having a higher mass. There may be a useful analogy to AGN in that their spectra are 
known to remain roughly constant on timescales of at least a few years (Markowitz et al. 
2003) while transient StBH systems undergo strong energy spectral changes on timescales 
of the order of a few days (McClintock et al. 2006 and references therein). It could be that 
the persistence of NGC 5408 X-l in the approximately same spectral state is related to it 
being of intermediate mass between StBHs and AGN. Alternatively, perhaps the source is 
in an accretion state, like the so-called ultraluminous state of Gladstone et al. (2009), that 
is just not well sampled by the behavior of Galactic StBHs. However, either conclusion will 
remain tentative until NGC 5408 X-l is observed in a state with a lower value of the power- 
law photon index, i.e., a state similar to the low/hard state in StBHs. Such an observation 
would clearly be important in order to make a more definitive identification of its QPOs in 
the context of those observed in StBHs. 
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Assuming the constancy of the energy spectral features, i.e., disk flux, T, with the 
increasing QPO centroid frequency (see Figure 5) is analogous to the saturation observed in 
the timing - spectral correlations in StBHs, one can estimate the minimum mass of the black 
hole in NGC 5408 X-l. If we are indeed seeing only the saturated portion of the complete 
timing - spectral correlation curve, then the minimum QPO frequency observed in NGC 
5408 X-l serves as an upper limit on the so-called transition frequency, i.e., the frequency 
beyond which the timing - spectral correlations tend to saturate. In order to derive a mass 
estimate, this particular frequency can be scaled to a reference StBH of known mass with 
a measured transition frequency {u trans ). We can use the transition frequency as defined by 
Shaposhnikov & Titarchuk (2009) in their fits of QPO and spectral parameters. We find a 
sample of three StBHs with mass and transition frequency measurements in Shaposhnikov & 
Titarchuk (2009) that are suitable for this purpose. Further, it is known that the same source 
can exhibit different tracks (timing - spectral correlation curves), i.e., different transition 
frequencies, in various outburst episodes. Given that we are interested in a lower limit on 
the mass of the black hole in NGC 5408 X-l, we consider the lowest transition frequency 
for a given source. The minimum mass of the black hole in NGC 5408 X-l, can then be 
calculated as v trana xMa,ssR e f erence /( minimum QPO frequency). The mass estimates using 
the three reference sources, GRO 1655-40, XTE J1550-564 and GX 339-4 a, re reported in 
Table 6 (~ 1000 M©), and show substantial overlap with those reported by S09. 


7. Summary 

The ULX in the irregular dwarf galaxy, NGC 5408 has been suggested to harbor an 
intermediate-mass black hole (Mass## 100 - a fewxlOOO M©) (S09). This mass estimate 
was strictly based on the assumption that the mHz QPOs seen from this source are the 
type-C analogs of the low-frequency QPOs (0.1-15 Hz) from StBHs. Here we have presented 
results from new observations of NGC 5408 X-l in which we searched for timing - spectral 
correlations similar to those often exhibited by the type-C LFQPOs from StBHs. Our analy- 
sis of multi-epoch XMM-Newton data from NGC 5408 X-l reveals that certain characteristic 
features of the power spectra, especially the QPO centroid frequency, changed significantly. 
However, the energy spectrum has remained roughly constant. These results can be inter- 
preted in two ways. This could be due to complete independence of the timing properties 
on the energy spectra, unlike StBHs, in which case mass scalings derived from the QPOs 
are likely to be problematic; or there is in fact a correlation, but w r e are seeing only the sat- 
urated part of the correlation behavior (constancy of the F/disk flux with increasing QPO 
frequency). Such saturation is often seen in StBHs (e.g., Vignarca et al. 2003). Assuming 
we are seeing this saturated portion of the correlation curve, we estimate the low r er limit on 



the mass to be 800 Mq. At least one observation with the source in a low/hard-like state 
(F < 2.0) is necessary to resolve the issue of whether the timing and spectral properties are 
correlated as in StBHs or not. 
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Figure 1A: / Left Panel]: Good time intervals (>5 ks) from EPIC-pn, MOS1 & 2 onboard 
XMM-Newton. For a given instrument, at a given time (X-axis), a finite value on the 
y-axis implies the instrument was continuously active for at least 5 ks. Also shown is a 
combined background (pn+MOSl+MOS2) light curve with flaring evident in some cases. 
The observation year is indicated at the top of each plot (Shown here are 2006, 2008 & 
2010A data: see Table 1). Time zero is indicated in secs since 1998.0 TT. [Right Panel]: 
Filtered light curves (pn+MOSl+MOS2) accounting for both instrumental dead times and 
flaring. The power spectra were derived from these cleaned light curves. These light curves 
were extracted using photons in the energy range of 0.3-10.0 keV. 
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Figure IB: [Left Panel]: Good time intervals (>5 ks) from EPIC-pn, MOS1 &: 2 onboard 
XMM-Newton. For a given instrument, at a given time (X-axis), a finite value on the 
y-axis implies the instrument was continuously active for at least 5 ks. Also shown is a 
combined background (pn+MOSl+MOS2) light curve with flaring evident in some cases. 
The observation year is indicated at the top of each plot (Shown here are 2010B, 2011A 
& 2011B data: Refer to Table 1). Time zero is indicated in secs since 1998.0 TT. / Right 
Panel]: Filtered light curves (pn+MOSl+MOS2) accounting for both instrumental dead 
times and flaring. The power spectra were derived from these cleaned light curves. These 
light curves were extracted using photons in the energy range of 0.3-10.0 keV. 
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Figure 2A: Shown are the power spectra derived from the first three of the six XMM-Newton 
observations (2006, 2008, 2010A: Table 1). For a given power spectrum the strength of a QPO 
is dependent on the energy band under consideration. [Left Panel]: The power spectrum using an 
energy band (shown at the top of each plot) in which the QPO is clearly detected. Also shown 
are the sizes of the typical error bars and the residuals (Data-Model) offset to an arbitrary value 
in each case. / Right Panel]: Same as the plots on the left panel. However, here we choose a fixed 
energy range of 1.0-10.0 keV. A consistent energy band across all the power spectra is required 
to unambiguously assess the behavior (especially, rms amplitude versus QPO centroid frequency: 
Fig.3) of these QPOs. 
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Figure 2B: Shown are the power spectra derived from the last three of the six XMM-Newton 
observations (2010B, 2011A, 2011B: Table 1). For a given power spectrum, the strength of a 
QPO is dependent on the energy band under consideration. / Left Panel]: The power spectrum 
using an energy band (shown at the top of each plot) in which the QPO is clearly detected. Also 
shown are the sizes of the typical error bars and the residuals (Data-Model) offset to an arbitrary 
value in each case. [ Right Panel]: Same as the plots on the left panel. However, here we choose 
a fixed energy range of 1.0-10.0 keV. A consistent energy band across all the power spectra is 
required to unambiguously assess the behavior (especially, rms amplitude versus QPO centroid 
frequency: Fig.3) of these QPOs. 
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Figure 3: The rms amplitude of the mHz QPOs detected from NGC 5408 X-l (Y-axis) is 
plotted against the QPO centroid frequency (X-axis). To study the dependence of this plot 
on the bandpass, we consider two different energy bands. [Left Panel]: Centroid frequency 
versus QPO amplitude using photons in the energy range from 1.0-10.0 keV. / Right Panel]: 
Same as the left panel, but a different band-pass of 0.5-5. 0 keV was used. The qualitative 
nature of the behavior does not change with the energy band under consideration. However, 
the rms amplitude of a given QPO seems to increase with increasing energy (at least in the 
two energy bands considered), i.e. , higher rms amplitude in 1.0-10.0 keV compared to 0. 5-5.0 
keV. 
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Figure 4: [Top Panel]: The X-ray energy spectrum (0.3-10.0 keV) of NGC 5408 X-l derived 
from the 2010B observation (Table 1). The data from all the three detectors is shown. Black 
corresponds to EPIC-pn, while green and red represent data from the MOSl and MOS2 
detectors, respectively. The spectrum was re-binned for clarity to have a significant detection 
of at least lOcr. However, no more than 10 neighboring bins were combined to achieve this 
level of significance, i.e,, setplot rebin 10 10 in XSPEC. [Middle Panel I: The ratio of the data 
to the model defined by tbabs* ( diskpn + apec + powerlaiu) is shown. [Bottom Panel;: The 
ratio of the data to the model defined by tabs* (bmc+apec) is shown. Both the ratios were 
re-binned to have a detection significance of at least 20 with no more than 20 neighboring 
bins combined, i.e., setplot rebin 20 20 in XSPEC. 
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Figure 5A: Timing-Spectral Correlations using the phenomenological model: 
diskpn + power law - Six observations/data points were used to search for timing-spectral 
correlations. [Top Left Panel]: The photon index of the power-law (Y-axis) is plotted against 
the centroid frequency of the QPO (X-axis). [Top Right Panel]: The disk contribution to 
the total flux (Y-axis) in the energy range of 0.3-10.0 keV is plotted against the centroid 
frequency of the QPO. [Bottom Left Panel]: Flux from the power-law component (Y-axis) is 
plotted against the QPO centroid frequency (X-axis). [Bottom Right. Panel]: For the sake of 
completeness, the disk temperature (Y-axis) is plotted against the QPO centroid frequency 
(X-axis). Within the present scheme of X-ray states in StBHs. these results are qualitatively 
consistent, with XGC 5408 X-l being in the steep power law state exhibiting type-C QPOs, 
with the characteristic timescales scaled down by a factor of ~ a fewxlO. 
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Figure 5B: Timing-Spectral Correlations using the BMC model - [Top Panels ]: 
The behavior of a typical stellar-mass black hole (1998 outburst of StBH XTE J1550-564) is 
shown (we plot data from Shaposhnikov & Titarchuk 2009). For clarity, the X-ray state of 
the source is highlighted using the following color scheme. Blue squares indicate the low/hard 
state, magenta corresponds to a state intermediate between the low/hard and the SPL state 
while the green triangles highlight the SPL state. From left to right, the dependence of 
the centroid frequency of the QPO, BMC normalization (disk flux) and the Comptonized 
fraction on the power law index are shown. The error bars are not indicated when their size 
is smaller than the size of the data point. [ Bottom Panel]: Similar plots using the data from 
NGC 5408 X-l. The results from the ULX NGC 5408 X-l are qualitatively consistent with 
the source being in an X-ray state very similar to the steep power-law state in stellar-mass 
black holes (Compare with the green triangles on the top panel). 









Table 1: Summary of the XMM-Newton observations of NGC 5408 X-l 


ObsID 

Tabs 1 (ks) 

Net Countrate 2 (cts s 1 ) 

Effective Exposure 3 (ks) 

0302900101 (2006) 

132.25 

1.26 ±0.04 

99.94 

0500750101 (2008) 

115.69 

1.19 ±0.04 

48.55 

0653380201 (2010A) 

128.91 

1.46 ±0.04 

80.68 

0653380301 (2010B) 

130.88 

1.40 ±0.04 

110.0 

0653380401 (2011A) 

121.02 

1.34 ±0.07 

107.75 

0653380501 (201 IB) 

126.37 

1.31 ±0.07 

98.66 


1 The total observation time. 

2 The average pn+MOS countrate. 

3 After accounting for flaring background and good time intervals. See §2 and §3 for details 
on filtering. 



30 


Table 2: Summary of the best-fitting model parameters for the power density spectra (PDS) 
of the ULX NGC 5408 X-l. Here, we use an energy band in which the QPO is detected with 
very high confidence. 


Dataset 

2006 

2008 

2010A 

2010B 

2011A 

201 IB 

Exposure 3, (ks) 

27.6x3 

15.5x2 

23.2x3 

21.0x4 

67.4x1 

55.0x1 

Energy Range b (keV) 

0.30-2.0 

1.0-10.0 

0.85-10.0 

1.0-4. 0 

1.10-8.0 

1.20-2.0 

Countrate 0 (cts s -1 ) 

1.17 ±0.03 

0.43 ±0.02 

0.71 ± 0.02 

0.50 ± 0.01 

0.41 ±0.02 

0.23 ±0.01 

c* 

1.95 ±0.01 

1.81 ±0.02 

1.95 ± 0.02 

2.00 ±0.01 

1.95 ±0.02 

1.92 ±0.02 

A* 

0.86 ± 1.05 

3.39 ±2.90 

9.14 ± 17.90 

(0.80 ± 2.30) X 10-° 4 

3.49 ± 10.61 

(2.11 ±3.47) X 10-° 3 

^Low* 

0.20 ±0.19 

0.09 ±0.13 

-0.25 ±0.31 

2.68 ± 0.98 

-0.04 ±0.42 

1.45 ±0.52 

v bend* (mHz) 

6.33 ±0.75 

6.34 ± 0.49 

10.26 ± 6.03 

30.63 ±20.20 

5.80 ±4.84 

7.15 ±8.67 

r High* 

6.74 ±3.34 

10.25 ±5.60 

1.60 ± 0.50 

0.25 ±0.17 

1.98 ± 1.83 

-0.09 ± 0.43 

Nqpo.P 

1.88 ±0.16 

3.51 ± 1.69 

0.63 ±0.15 

0.72 ±0.21 

2.10 ±0.56 

1.04 ±0.29 

i/o,i f (mHz) 

17.68 ± 0.77 

10.28 ±0.24 

40.40 ± 2.93 

38.81 ±2.01 

18.67 ±2.20 

19.43 ± 0.90 

Ais^ (mHz) 

15.48 ±2.05 

1.28 ±0.97 

28.06 ± 11.96 

18.75 ± 8.98 

21.98 ±6.23 

7.08 ± 3.60 

Nqpo,2 t 

- 

1.91 ±0.35 

- 

- 

- 

- 

1*0,2^ (mHz) 

- 

15.05 ± 1.68 

- 

- 

- 

- 

Ai/p (mHz) 

- 

11.10 ±3.05 

- 

- 

- 

- 

,\ 2 /dof 

(continuum 11 ) 

309.83/208 

(618.71/211) 

753.75/764 

(868.30/770) 

178.14/137 

(257.33/140) 

128.05/97 

(207.54/100) 

327.11/329 

(669.58/332) 

268.86/267 

(543.03/270) 


a 

The good time intervals were broken into smaller intervals to improve the signai-to-noise in the power density spectra (PDS). 
For a given observation, the size of each segment x number of such segments is shown. 

b The power spectrum was derived using all the photons in this energy range. For a given observation, this is 
an energy range in which the QPO was detected with a high significance. 

c The countrate in the bandpass shown in the second row. 

* We fit the continuum with a bending power-law model described as follows: 


Continuum — C + 


A.U 1 Low 


High 1 Low 


where, V Low and T High are the low and high frequency slopes, respectively, and Vbcnd iS the bend frequency, 
t We model the QPOs with a Lorentzian. The functional form is as follows: 


QPO = 


i+ ( 


Nqpq 


T 


where, uq Is the centroid frequency and Ai/q is the FWHM of the QPO feature. 
d The x 2 /dof for the continuum are shown in braces. 



-31 - 


Table 3: Summary of the best-fitting model parameters for the power density spectra (PDS) 
(using only photons in the energy range of 1.0-10.0 keV) of the ULX NGC 5408 X-l 


Dataset 

2006 

2008 

2010A 

2010B 

2011A 

201 IB 

Exposure 4 * (ks) 

27.6x3 

15.5x2 

23.2x3 

21.0x4 

67.4x1 

55.0x1 

Energy Range (keV) 

1.0-10.0 

1.0-10,0 

1.0-10.0 

1.0-10.0 

1.0-10.0 

1.0-10.0 

Countrate l> (cts s -1 ) 

0.43 ±0.01 

0.43 ± 0.02 

0.55 ± 0.02 

0.53 ±0.01 

0.49 ± 0.02 

0.49 ± 0.02 

c* 

2.00 ±0.01 

1.81 ±0.02 

1.90 ±0.04 

1.99 ±0.01 

1.94 ±0.01 

1.92 ±0.01 

A* 

0.80 ± 0.89 

3.39 ± 2.90 

0.01 ±0.01 

(1.48 ±3.88) X 10-° 4 

4.16 ± 11.49 

2.96 ± 9.50 

F Low * 

0.20 ±0.17 

0.09 ±0.13 

1.12 ±0.30 

2.65 ± 1.06 

-0.072 ± 0.39 

-0.05 ± 0.56 

Vbend* (mHz) 

7.89 ± 0.73 

6.34 ± 0.49 

7.12 ±6.13 

36.39 ±21.26 

4.96 ± 1.48 

8.30 ± 1.15 

T High * 

10.02 ± 11.47 

10.25 ±5.60 

-0.41 ±0.46 

0.25 ±0.14 

4.42 ± 7.78 

9.87 ± 12.11 

Ngpo.d 

2.49 ±0.18 

3.51 ± 1.69 

0.64 ±0.15 

0.67 ±0.21 

2.40 ±0.19 

2.64 ± 0.27 

i/o.i* (mHz) 

16.81 ± 1.05 

10.28 ±0.24 

39.07 ± 1.70 

39.07 ± 1.82 

18.49 ± 1.28 

15.04 ±2.04 

AiT (mHz) 

19.51 ±2.16 

1.28 ±0.97 

15.15 ±6.44 

16.04 ±8.89 

21.86 ±2.60 

22.83 ± 2.94 

RMSgpo.l^ 

42. 12 ±2.83 

12.81 ±5.75 

16.64 ± 4.05 

17.84 ±5.68 

41.00 ±3.05 

43.94 ± 3.72 

NQPO,2 t 

- 

1.91 ±0.35 

- 

- 

1.31 ±1.02 

- 

i^ 0 , 2 ^ (mHz) 

- 

15.05 ± 1.68 

- 

- 

8.40 ±1.22 

- 

Aiqt (mHz) 

- 

11.10 ±3.05 

- 

- 

- 

- 

RMSqpo, 2* 

- 

27.82 ± 4.64 

- 

- 

3.09 ±4.32 

- 

X 2 /dof 

(continuum** ) 

268.34/208 

(755.28/211) 

753.75/764 

(868.30/770) 

200.60/137 

(242.10/140) 

135.27/97 

(190.64/100) 

350.16/326 

(851.45/332) 

291.45/266 

(677.00/269) 


The good time intervals were broken into smaller intervals to improve the signal- to- noise in the power density spectra (PDS). 
For a given observation, the size of each segment x number of such segments is shown. Same time intervals as in Table 1 were 
used here. 

b To systematically study the QPO properties, we reduced the power spectra in a consistent energy band of 
1.0-10.0 keV. 

c The count rate in the bandpass of 1.0-10.0 keV. 

* We fit the continuum with a bending power-law model described as follows: 


Continuum — C + 


1 + 


Au ^ Loiu 

f » 

v u be.nd ' 


where, Vi ow and F High are the low and high frequency slopes, respectively, and Vbend is the bend frequency, 
t We model the QPOs with a Lorentzian. The functional form is as follows: 


QPO = — i RMSqpo (integrated f rom — oc to -f oo) — x 100 


where, vq is the centroid frequency and A*/o is the FWHM of the QPO feature. 
The ,y 2 /dof for the continuum are shown in braces. 
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Table 4: Summary of the energy spectral modeling of NGC 5408 X-l. [Top Panel]: Best- 
fitting parameters for the tbabs* (diskpn+apec+pow) model. [Bottom Panel]: Best-fitting 
parameters for the tbabs* (bmc+apec) model. 


tbabs * ( diskpn+apec+pow P 


Dataset 

ll H 

'-pb 

1 max 

kT c 

plasma 

r d 

Flux e 0.3— 10.0/ee V 

Flux f Diak 

X 2 /dof 

2006 

13.84+g-™ 

n 144~*~ 9005 

U ‘ A44 -0.005 

0- 90 to.06 

2 64+°*°4 

Z ‘°4-0.04 

3.92l21 3 o x 10-12 

1.72i°4|xl°- 12 

1035.02/837 

2008 

14.35ij.43 

0 i37+ () 01 ° 
u,AO ' -0.009 

0 - 8 ito.S! 

o kc+0.06 
Z ‘ OO -0.06 

3.73lgJ?xlO- 12 

1.54l°;|‘x io-» 

614.50.02/603 

2010A 

1 ft '7'7~t-0.86 
Aa *» * -0.82 

n i^o +() 007 
U,AOU -0 .006 

0.921°;°? 

o kc+ 0.04 
z *°°-0.04 

4.36l°;‘|xi0- 12 

1.65lg;igxl0-» 2 

1072.47/888 

2010B 

12.48l°;?| 

0 160+ 0 007 
U.1DU_ () 007 

o.86l°;°| 

O kc+ 0.04 
^* Ot> -G.04 

3.90l“;»xl0~ 12 

i.27ig;iSxio-« 

1110.21/929 

2011 A 

14.40i°;“ 

n i4s+ 0007 

U- 140-0 .006 

0-84l°;“ 

o ci +0.04 
‘ 4 * OA -0.04 

4-15lSilxlO- 12 

1.62l“;^xl0- 12 

915.33/838 

2011B 

13.051°;“ 

n i C7+0.008 
u,AO/ -0.008 

0.89l°;J° 

O co+0.05 
4.0 ^_o > o5 

3.73l5j;yxl°- 12 

!.39l°;“xl0- 12 

855.71/779 


tabs * ( bmc+apec ft 


Dataset 


kTK diak 

a=r h -1 

f 

BMC 

kT c 

plasma 

X 2 /dof 

2006 

Q 17+0.95 
y ‘ A * —0.80 

nil q+o. 003 
U,li °-0 .003 

1 fifi+0-03 
ado -0.03 

o-ssio:” 

3.03l°;?? xl0“ 5 

0 Q4 - * -0 05 
Uy4 -0.06 

947.01/837 

2008 

10.46li;|| 

A 11 9 +O.OO 6 
U11Z -0.006 

1 K 7 +O.O 6 
1,0 * -0.06 

0-39l°;°5 

2.98-1“;“ X10- 5 

H qo+9.07 

U.04 — 005 

598.82/603 

2010A 

Q 11 +105 
11 -0.95 

0 11Q+0.00 4 

U.iiy _ 0 004 

1 ioft+003 

A 00 — 0-03 

0 4 *;+ 0-02 
U,4 ' :) -0.03 

3.171°;!* xio~ 5 

0 qk+005 

u-yo_o.o5 

1022.89/888 

2010B 

'•4o_.o-.38 

A 1 9O+0.004 
U.iZ4 _ 0 o 0 4 

1 cq+0.03 
A - oy -0.03 

n ko+ 0.02 

.03 

2.65l°;?fxl0 -5 

0 Q 1 +9 - 06 

UyA -0.05 

1041.68/929 

2011A 

Q 07 + 1.12 
y,,5/ -1.06 

0 1 20 +() 004 
U. A zU_o o 04 

1 «o+0.04 
AD,3 -0.04 

A 4K+0.03 
°‘ 45 -0.03 

3.0°i°;|xl0- 5 

a q-7+0.06 

°- 87 -0.04 

845.07.07/838 

2011B 

o 41 + A-2A 
8 ’ 4A — 1.02 

A 1 9O+0.005 

U. xZo_ q Q05 

1 ck+0.04 
a ' oo -0.04 

u * iD -0.03 

2.67l°;ixl0- 6 

0 Q4+0-08 
U ‘ y4 -0.08 

829.33/779 


Total column density of hydrogen along the line of sight including the Galactic extinction (in units of 10 2o cm~ 2 ). We used 
the tbabs model in XSPEC. b Accretion disk temperature in keV. We used the diskpn model in XSPEC. The inner radius of 
the disk was fixed at 6GM/c 2 . c The temperature of the surrounding plasma in keV. We used the apec model in XSPEC. The 
abundances were fixed at the solar value. 4 The photon index of the power law. e The total unabsorbed X-ray flux (in units 
of ergs cm" 2 s" 1 ) in the energy range of 0.3-10.0 keV. * The disk contribution to the total X-ray flux (in units of ergs cm -2 
s” 1 ) in the energy range of 0.3-10.0 keV. 

g The color temperature of the disk blackbody spectrum (in keV). h Spectral index of the power law portion of the energy 
spectrum. s The Comptonized fraction: fraction of the input blackbody photons that are Compfconized by the bulk motion of 
the in-falling material. } The normalization of the disk blackbody spectrum (in units of (L/10 39 erg s"" i )/(d/10 kpc) 2 ). This 
is an indicator of the disk flux. 
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Table 5: Comparing the PDS and the energy spectral properties of NGC 5408 X-l with the 
steep power-law (SPL) state in stellar-mass black holes 


X-ray State 

Definition 8, 

Steep Power Law State 
(SPL) 

Presence of power-law component with T > 2.4 
Power continuum: r 6 < 0.15 

Either disk fraction e < 80% and 0.1-30 Hz QPOs present with rms amplitude > 0.01 
or disk fraction 0 < 50% with no QPOs present 

X-ray State 

Data from NGC 5408 X-l d 

Steep Power Law State ? 
(SPL) 

Presence of power-law component with P in the range 2. 5-2. 7 
Power continuum: r* in the range « 0.36 - 0.43 
Disk fraction® in the range ~ 30-50% and 
10-40 mHz QPOs present with rms amplitude in the range 0.10-0.45 


a 

In case of stellar-mass black holes, the definition is using data in the energy range of 2-20 keV. 
b Total rms power integrated over 0.1-10.0 Hz. 
c Fraction of the total 2-20 keV unabsorbed flux. 

li In case of NGC 5408 X-l, we use the XMM-Newton bandpass of 0.3-10.0 keV. 

* We integrate the total power (continuum + QPOs) in the frequency range of 0.001 - 0.5 Hz (same as the 
frequency range used for power spectra fitting). The power spectra were derived using all the photons in 
the energy range of 1.0-10.0 keV (see Table 3 for the parameters). 

0 Fraction of the total 0.3-10 keV unabsorbed flux. 
f Fraction of the total 0.3-10 keV unabsorbed flux. 
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Table 6: Mass estimate of the black hole in NGC 5408 X-l 


Reference 

Mass(M 0 ) 

L'trans (Hz) 

Min. Mass^ GC . 5408 (MQ) 

References 

GRO J1 655-40 

6.3 ±0.5 

3.0 ±0.1 

1720 

1,4 

XTE J1550-564 

9.5 ±1.1 

1.84 ±0.07 

1490 

2,4 

GX 339-4 

> 6 

1.4 ±0.2 

820 

3,4 


a 


The transition frequencies were estimated by Shaposhnikov et al. (2009). The fitting 
functions for the timing-spectral correlations are described therein. 

b 

Minimum mass of the black hole in the ULX = (^ r aras x niass)fl,./erence5tB/f/10.28 mHz, 
where 10.28 mHz is the minimum QPO frequency detected from NGC 5408 X-l. We 
consider the error bars on the individual parameters and report the lower value of the mass. 
References - (1) Greene et al. 2001; (2) Orosz et al. 2002 (3) Munoz-Darias et al. 2008 (4) 
Shaposhnikov et al. 2009. 



